The Singularity Threshold of the Nonlinear Sigma Model Using 3D Adaptive Mesh 
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Numerical solutions to the nonlinear sigma model (NLSM), a wave map from 3 + 1 Minkowski 
space to S 3 , are computed in three spatial dimensions (3D) using adaptive mesh refinement (AMR). 
For initial data with compact support the model is known to have two regimes, one in which 
regular initial data forms a singularity and another in which the energy is dispersed to infinity. The 
transition between these regimes has been shown in spherical symmetry to demonstrate threshold 
behavior similar to that between black hole formation and dispersal in gravitating theories. Here, 
I generalize the result by removing the assumption of spherical symmetry. The evolutions suggest 
that the spherically symmetric critical solution remains an intermediate attractor separating the 
two end states. 



Nonlinear sigma models have attracted the attention 
both of physicists and mathematicians. For physicists, 
the models represent the possibility of describing par- 
ticles with a field theory, while mathematicians study 
singularity formation in geometrically motived nonlinear 
models. Recently, a nonlinear sigma model has attracted 
the attention of relativists because it demonstrates be- 
havior similar to black hole critical phenomena discov- 
ered by Choptuik Q. 

Studying the gravitational collapse of a spherically 
symmetric scalar field, he found that initial data with 
large energy formed a black hole while that for small 
initial energy dispersed its energy to infinity. By tun- 
ing the initial energy, he approached the so-called criti- 
cal regime which separates these two end states. In this 
regime, solutions approach a unique, universal solution 
which demonstrates self-similarity. 

The model studied here serves as a useful toy model 
for such threshold behavior because it also possesses two 
stable end-states, dispersal and singularity formation. 
These end-states serve as analogues for the gravitating 
case in which initial data can either disperse or form 
a black hole. Similar to Choptuik's gravitating model, 
spherically symmetric studies of this model have found a 
universal, self-similar critical solution (|, ||. The utility 
of the nonlinear sigma model arises from its relative sim- 
plicity compared to any of a variety of gravitating mod- 
els. Perhaps the study of such nonlinear, flat systems 
will guide the way to further insight in the gravitating 
models. 

While questions about the nature of black hole crit- 
ical behavior remain, many unexplored "experiments" 
have yet to be conducted. While many matter models 
have been studied (for a review see Q), generally these 
studies have been restricted to spherical symmetry. The 
only work to model fully nonlinear collapse in less sym- 
metry is that of Abrahams and Evans Q who studied 
vacuum gravitational collapse in axisymmetry. One rea- 



son for the scarcity of multi-dimensional evolutions of 
the critical regime is the extreme resolution demanded 
by approach to scale invariance. Perturbation methods 
are a complementary approach with recent work [tJ 
suggesting that spherical critical solutions in the scalar 
field and perfect fluid cases extend to the nonspherical 
regime. While work is underway to duplicate and extend 
the work of Abrahams and Evans in axisymmetry j|] , it 
makes sense to look ahead to a simple model in three 
spatial dimensions. 

To obtain the resolution required for evolving self- 
similar solutions, Choptuik developed a computational 
infrastructure in one spatial dimension which dynami- 
cally and locally adds numerical resolution where needed 
using adaptive mesh refinement (AMR) . With this infras- 
tructure, fine subgrids are added and subtracted to the 
computational domain providing resolution only where 
and when needed. Because the computational cost scales 
as a power law in the spatial dimension (i.e. doubling the 
resolution of a rf-dimensional evolution requires a factor 
2 d+1 more work), it is expected that AMR will be abso- 
lutely crucial in higher dimensions for well-resolved evo- 
lutions of interest, in particular black hole critical phe- 
nomena or black hole collisions. 



Constructing such a gravitating model with AMR is an 
ambitious project. Instead I report on the construction 
of a 3D AMR code of the nonlinear sigma model, which 
genuinely requires AMR and holds physical interest. 

The NLSM Model: The nonlinear sigma model 
studied here represents a mapping from the base space 
of 3 + 1 Minkowski to a target space of S 3 . In spherical 
symmetry, it is common to choose the hedgehog ansatz 
for the map reducing the dynamics to that of a single 
spherically symmetric field x(r, t). Here a simple gener- 
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alization of this ansatz is chosen 



( sin x(x, y,z,t) sin 8 sin ip \ 
sin x(a;, y,z,t) sin cos y> 
sinx(a;, y, 2, i) cos (9 
cosx(^,y,2,t) 



(1) 



where 9 and y are the usual spatial angles. The dynamics 
reduce to the scalar field x(a;,y, z, t) which satisfies the 
equation of motion 



X = X,: 



X,: 



X„ 



sin 2x 



(2) 



where commas indicate partial differentiation with re- 
spect to subscripted coordinates, an overdot denotes 
d/dt, and r = ^/ x 2 + y 2 + z 2 . The equation of mo- 
tion (0) implies the regularity condition x(0,0, 0,i) = 
which is enforced by the evolution procedure. This ansatz 
requires that the origin is singled out as a special point. 
More elegant generalizations may be considered in the 
future. 

The energy density of the map is given by 



1 

P= 2 



(x) 2 + (X,,) 2 + (X,y'f + (X, Z ) 2 



sin 2 x 



The angular momentum densities are given by [^| 
M'" = I d 3 x (T^x" - T°"x^) , 



(3) 



(4) 



so that the z-component of the angular momentum, for 
example, is 



J z = J d A x M xy = J d A x X (yx,x ~ *X,v) ■ (5) 

While the map allows for the possibility of a texture 
charge associated with the third homotopy group, only 
initial data with zero charge is considered here. The 
model requires initial data xi x i 2/> z > 0) an d xOej Z/> 0) 
be specified at the initial time. 

Various types of initial data have been implemented 
and are described in Table |[ Some of these families are 
defined in terms of a generalized Gaussian pulse defined 
by 



G(x,y,z)=Ae^-V 2 ^, 
where f is a generalized radial coordinate 

f = \fe x (x 



(6) 



x c f + £ y (y- y c f + (z- Z c f. (7) 



Such a pulse depends on parameters: amplitude A, shell 
radius R, pulse width S, pulse center (x Cl y c , z c ), and 
skewing factors e x and e y . For e x ^ 1 =/= e y such a 
pulse has elliptic cross section. Family (a) represents 
a single pulse for which the parameter v takes the val- 
ues {—1,0,-1-1} for an approximately out-going, time- 
symmetric, or approximately in-going pulse. The angular 





Description 




X(x,y,z,0) 


a 


Ellipsoid 


G 


or 

+n z (yG, x -xG, y ) 


b 


Two pulses 


G\ + G2 




c 


Toroid 


Ae -z 2 /S 2 e -{e** 2 +e V V 2 y/ S2 






TABLE I: List of various initial data families. For fami- 
lies (a)-(c) both the field x( x i Vi z i 0) and its time derivative 
x(x,y,z,0) are shown in terms of various parameters. The 
terms G, Gi, and G2 represent unique Gaussian pulses as de- 
fined in Eq. (^|). In family (b), the parameters vi and V2 are 
the respective velocities of the two pulses, generally chosen to 
have a grazing collision. 



momentum of the pulse about the z-axis is proportional 
to the parameter Sl z as well as to (e x — e y ) 2 . 

Numerical method: RNPL j[o) is used to develop 
and debug a stable and convergent unigrid code which 
solves Eq. (||) using finite differences and iterative Crank- 
Nicholson. The RNPL generated update procedure is 
called from the AMR code to evolve any given fine grid. 
The AMR implementation follows that of Berger and 
Oliger with some simplifications: 

(1) Instead of estimating truncation error, I use the 
energy density as a criterion of refinement. In particular, 
normalizing the energy density p by the grid resolution 
1/h, the refinement criterion is h 2 p > e, where e is a user- 
specified threshold. This simplification is less general 
than truncation error but provides an easily computable 
and smooth function from which to estimate resolution 
requirements. 

(2) Fine grids are completely contained within their 
parents. Grids do not overlap or abut other grids at 
the same level. This restriction reduces the maximum 
obtainable efficiency but provides for considerable sim- 
plification. 

(3) Fine grids are created strictly aligned with parent 
grids with no rotation. The transformations necessary to 
relax this would be especially onerous in the gravitating 
case. 

(4) The ratio of refinement between parent and child 
grids is constrained to be an even integer. 

The first two of these restrictions appear especially 
suited for the case of central collapse studied here, and 
the intention is to relax them for more general problems. 

The usefulness of AMR depends obviously on whether 
it produces correct solutions and whether it allows for 
high resolution with only proportional work. Starting 
with the unigrid code, convergence and energy conser- 
vation were confirmed. The AMR results were then 
checked against high resolution unigrid results. Perhaps 
a stronger test is that this code finds the same results 
for spherically symmetric initial data as that in (as 
demonstrated in Fig. Fy). As an example of the benefit 
of AMR consider Fig. g. A rough estimate of the com- 




FIG. 1: Approach of near-critical evolution to self similarity. 
The numerical evolution is shown (circles) for spherically sym- 
metric initial data at times near the collapse time. The initial 
data is a member of family (a) in Table | with e x = e y = 1 
and Q. z = 0. The data represents a (x > 0, y — 0, z = 0) cut 
with every point shown. The excited n = 1 self-similar solu- 
tion is shown (solid line) with r = In \T* — T\ where T* is the 
time of collapse (so that collapse occurs at r — > -co). The 
collapse time of the n = 1 solution is chosen so that the two 
solutions coincide for the first frame only. That they coincide 
for the other frames indicates the approach to the self-similar 
solution. 



putational work for the entire evolution is 2 8 , in units 
where the work for a unigrid evolution at the coarsest 
resolution is unity. In comparison, to achieve uniform 
resolution equal to that of the finest sub-grid would re- 
quire work equal to (2 9 ) 3+1 , a factor 2 26 more work than 
that with AMR. 

Results: As discussed in Q, solutions generically 
evolve to one of two stable end states. For large initial 
data, the energy density collapses toward the origin sug- 
gesting the formation of a singularity. As in the strictly 
spherical case, the usual caveat applies that we only have 
indications, not proof, that a singularity is forming. For 
small initial data, the energy density implodes through 
the origin and then disperses to infinity. Removing the 
restriction to spherical symmetry has not revealed any 
additional stable end states. 

The aim is to investigate the region of solution space 
separating these two end states. In the jargon of non- 
linear dynamics, we look at the so-called critical surface 
occurring between the two basins of attraction. To do 
so, one chooses a 1-parameter family of initial data with 
the property that data with small parameter disperses 
while large parameter initial data forms a singularity. By 
tuning the free parameter according to the resulting end 
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FIG. 2: Snapshots of slightly sub-critical evolution. The 
field x(t> x >V> z = 0) is shown at four times for a family of 
initial data with angular momentum about the z-axis. For 
clarity, only every other point in the x and y directions is 
shown (every fourth point). At the final time, all ten grids 
are shown with a refinement factor of 2. Greater detail of this 
last frame is shown in Fig. ^j. The initial data is a member of 
family (a) from Table | with parameters A w 1.3359, R = 8, 
8 = 3, x c = y c = z c = 0, e x = 0.5, e v = 1, u = 0, and 
Q z = 0.4. The dimensionless ratio of the angular momentum 
to the energy squared is J/E 2 — 0.0025. 




FIG. 3: Detail of the peak region from the final frame of 
Fig. [2|. Only the finest four grids are shown (every point is 
shown) along with the coarsest grid in order to compare the 
resolution achieved with respect to the coarse grid. The finest 
grid has a resolution a factor of 2 9 finer than the coarse grid. 
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state of the evolution, the critical surface is approached. 

The previous results in spherical symmetry indicated 
that all initial data of compact support, when tuned, 
approached the same universal, self-similar solution [Q. 
This critical solution happens to be one of a family of so- 
lutions found by Aminneborg and Bergstrom by assum- 
ing self-similarity and solving the resulting ODE jL2| . A 
linear perturbation analysis revealed the n = 1 solution 
to have the requisite single unstable mode necessary to 
be a proper critical solution (a so-called intermediate at- 
tractor) . 

This solution is shown in coordinates adapted to its 
self-similarity in Fig. Given a self-similar solution 
which takes the form f(r,t) = f(r/\T* —T\), one can 
recast it as g (bar — In \T* — T\). This form makes ap- 
parent that the solution executes linear motion in log 
space and log time, where one is free to set the collapse 
time T*. Thus, in Fig. [j] the self-similar solution travels 
leftward retaining its shape. 

Along with the n = 1 solution is shown a cut along the 
x-axis of a near-critical evolution from spherically sym- 
metric initial data. Such an evolution would be expected 
to agree with previous results in spherical symmetry, as it 
appears to do in the plot. The collapse time of the n = 1 
solution was chosen so that the two solutions agree in 
the first frame. The n — 1 solution is then completely 
determined at other times by translation in lnr. 

It was also observed in 0] that whenever the range 
of x exceeded it at a given time, a singularity would 
eventually form in its future. However, evolutions of 
non-spherically symmetric families do not obey such a 
condition, having exceeded ir while remaining nonsingu- 
lar. One still finds solutions which disperse {sub- critical) 
separated from those which form a singularity {super- 
critical) . 

For such non-symmetric initial data, the question is 
whether some different critical solution appears. If the 
n = 1 solution has non-symmetric unstable modes not 
present in spherical symmetry, then a change in criti- 
cal solution would be expected. However, all such non- 
symmetric variations of the types of initial data in Table || 
have failed to find a different critical solution. 

An example from initial data with angular momentum 
is shown in Fig. || The pulse evolves in a nontrivial man- 
ner eventually dispersing much of its energy toward the 
edges of the grid. However, near the origin, a spherically 
symmetric waveform emerges which propagates toward 
the origin. In this case, the waveform has negative am- 
plitude (the model is invariant with respect to the trans- 
formation x ~X)- Further detail of this waveform is 
shown in Fig. ^. The upward sweep at the origin corre- 
sponds to the regularity constraint mentioned previously 
that x(0,0,0,f) = 0. 

Another example of initial data with angular momen- 
tum is shown in Fig. |^. Two pulses are sent towards one 
another and, for sufficiently large initial amplitude, sin- 



m 

® 


,7 = 0.2 


(7 = 4 f 

/ ^ 

7 


cr = 8 r 

i i i 1 i i i 



LC 



X 



LO 



-4 -2 
In r 

FIG. 4: Demonstration of approach to a spherically sym- 
metric solution at criticality. Initially asymmetric initial data 
is tuned and slightly sub-critical solutions are shown above 
for different families. Here a is a measure of the asymmetry 
(<r = corresponds to spherical symmetry) which partially de- 
fines family (a) in Table f| via e x = 1 + c and t v = 1 — a. For 
each value of a, three slices are shown: (x > 0, 0, 0) (circles), 
(0, y > 0,0) (triangles), and (0, 0, z > 0) (crosses). Their 
general agreement near the origin demonstrates approach to 
spherical symmetry. Also shown (line) is the n = 1 self-similar 
solution (with collapse time adjusted to match the solution). 



gularity formation occurs. Four frames of a near critical 
evolution show the region near the origin approaching 
the spherically symmetric, self-similar solution. 

A more systematic exploration away from spherical 
symmetry is shown in Fig. ||. Here, a parameterizes the 
asymmetry of the initial data, and near-critical solutions 
are shown for four different values of a. Furthermore, 
for each case the region near the origin (the region which 
approaches self-similarity) is shown along the three dif- 
ferent axes. Their agreement suggests that this region 
approaches spherical symmetry. 

There remains, however, the possibility that some 
other spherically symmetric self-similar solution becomes 
the critical solution away from spherical symmetry. 
Though there is no evidence of any transition, the various 



excited self-similar solutions of 1 12 are very similar near 
the origin consisting of a peak followed by oscillations 
about x — 7r /2- To distinguish among them, a pertur- 
bation analysis allowing for non-symmetric modes would 
be appropriate, with the critical solution being the one 
with a single unstable mode. 

The evidence therefore suggests that the n = 1 self- 
similar solution remains the intermediate attractor away 
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FIG. 5: Snapshots of a slightly sub-critical collision of two 
grazing pulses. The field \ (t, x iV, z = 0) is shown in the x — 
y plane at four times for two pulses with linear momentum 
in ±x directions. At late times, the region near the origin 
(the fine grids in the center of the final two frames) rapidly 
approaches spherical symmetry and self-similarity. The initial 
data consists of family (b) from Table [i] with parameters A\ — 
A 2 = 11.5, Ri = R 2 = 0, 5i = <5 2 = 0.5, (x cl , y cl , z c i) = 

{ — Xc2, ~yc2, Zci) — (1.4, 0.7, 0), txl = €x2 = 1, €yl = <Ly2 = 1, 

and v\ — —V2 = 1. The dimensionless ratio of the angular 
momentum to the energy squared is J/E 2 — 0.0008. 



from spherical symmetry, but a perturbation analysis 
should settle the issue. Further work will likely be di- 
rected towards (1) using the AMR infrastructure with a 
gravitating model, (2) relaxing the noted simplifications 
in the clustering, and (3) distributing the grids using MPI 
along the lines of [ p^[ . 
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